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In the face of uncertain biological response to climate change and the many critiques concern¬ 
ing model complexity it is increasingly important to develop predictive mechanistic frame¬ 
works that capture the dominant features of ecological communities and their dependencies 
on environmental factors. This is particularly important for critical global processes such as 
biomass changes, carbon export, and biogenic climate feedback Past efforts have success¬ 
fully understood a broad spectrum of plant and community traits across a range of biological 
diversity and body size, including tree size distributions and maximum tree height, from me¬ 
chanical, hydrodynamic, and resource constraints Recently it was shown that global 
scaling relationships for net primary productivity are correlated with local meteorology and 
the overall biomass density within a forest Along with previous efforts, this highlights 
the connection between widely observed allometric relationships and predictive ecology. An 
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emerging goal of ecological theory is to gain maximum predictive power with the least num¬ 
ber of parameters Here we show that the explicit dependence of such critical quantities 
can he systematically predicted knowing just the size of the largest tree. This is supported 
hy data showing that forests converge to our predictions as they mature. Since maximum 
tree size can he calculated from local meteorology this provides a general framework for 
predicting the generic structure of forests from local environmental parameters thereby ad¬ 
dressing a range of critical Earth-system questions. 

Terrestrial net primary produetivity (NPP) is one of the key proeesses affeeting the overall 
elimate system, and given its dependenee on loeal meteorology, defines one of the most important 
global feedbaek processes Classically it has been shown that terrestrial biomass production 

correlates with precipitation and temperature but more recent efforts have pointed to existing 
biomass density as the strongest determinant of NPP and highlight the indirect consequences of 
the local environment potentially complicating the prediction of the biological response to climate 
change Here we show that biomass density and net primary productivity can both be determined 
from maximum tree size which in turn has been shown to be a calculable complex function of local 
climate and consequently a strong predictor of biomass density, as supported by empirical evi¬ 
dence Through these relationships we also predict the observed relationship between NPP and 
biomass density, where these two features taken together represent the dominant climate feedbacks 
either through carbon uptake or the biological effects on atmospheric water and heat fluxes 

In our derivations we first express all scaling relationships in terms of general exponents 
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without specifying their value. Although many of these can themselves be derived from an under¬ 
lying theoretical framework, we use both theoretical as well as empirical values for these exponents 
in making predictions. These are compared with data spanning the entire continental United States 
representing a large number of whole forests subdivided into over 9x10^ one km^ size grid cells 
covering all 36 ecoregions We discuss how these predictions depend on differences in scaling 
relationships for a given species or environment and show how forest maturity affects the observed 
scaling results. 

Recent work has shown that the general spacing of trees can be derived by combining mor¬ 
tality processes and the geometric scaling of vascular plants with the idea that the resource use per 
unit area is approximately constant in a given forest Consequently, the average distance, dk, 
between individual trees within a given size class is predicted to scale with trunk radius: 

4 = ( 1 ) 

where k denotes a linearly binned size class of trees^^’ ci is a constant, and is the trunk radius, 
which is related to tree height, hk, and the mass of the tree, nik, by where 

Ch and Cm are also constants, and au and am are scaling exponents In this relationship ad is 
the scaling exponent which, from theoretical and empirical work, is expected to be close to 1 
Its possible variation could depend on a variety of factors from the scaling of overall metabolic 
rate to that of the canopy. Since equation 1 is based on the idea that resource consumption per 
unit area is constant this further implies that ci should depend on resource availability. This can be 
shown by noting that the largest tree in the forest completely dominates its area and does not share 
resources with any neighboring tree (if it did then it would be possible for a larger tree to exist) 
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In a similar fashion, the radial extent of the roots of the largest tree, rroot,max, must also define the 
separation of trees, in which case we have for the largest tree 

d-max — ^root,max max 

leading to Ci = . In previous work it was shown that the root radius is related to the trunk 

'^max 

radius, = Croot'^max^ which in turn can be determined by the local resource environment 
Previous work has shown that the scaling is given by = 2/3 and c^ooi = 

where the quantities /Ss, 1^, r^, and n are all basic normalization constants of scaling whose 
precise definition can be found in Ref. Through its dependence on r^ax , ci is a complicated 
but calculable function of solar radiation, temperature, humidity, precipitation, and altitude. We 
can test this result using the distributions in Ref. where the largest tree has a trunk radius of 
54.5 cm. Using = 2/3, ad = I and the above expression for Croot, we find ci = 73, in 
excellent agreement with the corresponding published best fit range of Ci = 63 to 78 from Ref. 
These results illustrate that the absolute distance between trees within any given size class can be 
determined solely from the size of the largest tree. 

Previous work has also shown that the number of trees of a given size scales with trunk 
radius as 

nk = c^r‘^"- (3) 

where, by definition, there is only a single largest tree (so Umax = 1 ); C 5 is a constant given by 
C 5 = with On = — 2 from previous work. This relationship explicitly shows that the overall 
properties of the entire forest are determined by the size of the largest individual. Because these re- 
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suits define the overall distribution of trees we can further show that numerous other features of the 
forest can be straightforwardly determined from the maximum size, notably the average biomass 
density and net primary productivity, both of which are of critical importance to atmospheric feed¬ 
back processes and climate dynamics. We first show how whole forest totals depend on maximum 
size and from these we can calculate various macro-scale averages of the forest. 


The mass of an individual tree within a given size class is given by ruk = so the 

total mass of all trees in the forest is 


1+Q;n + 1/ Oi-n 


Mtnt — 


UkTUkdr = 


max 


r, 


l+Cln + l/ Oin 
0 


'ro 


l/a„ 


1 -f -f \joL^ 


(4) 


^1 + 1/ OLn 


which means that for large maximum size the total forest mass scales as Mtot oc Tmax assuming 
that 1 + an + l/cim > 0 as supported by the existing exponent values discussed later. Given that 
the metabolic rate of an individual tree is given by the total metabolic rate of the 

entire forest is 


Pf'max / ^l+an+OlB _ \ 

Btot = / nMr = 5or-- ^ ° - (5) 

Jro \ 1+an + aB J 

leading to Btot oc . Taken together with equation 4 this implies that whole forest metabolism 
scales with whole forest mass with an exponent of ■ A common expectation for the scaling 

of metabolism with step radius is q;b = 2 (equivalent to scaling as and it is also 

commonly assumed that = 3/8 in which case whole forest metabolism scales with whole 
forest mass following a power of 9/11 which is larger than the single tree scaling of 3/4, but still 
less than linear. 


The total area required to observe a full distribution of tree sizes can be found by first noting 
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that the largest tree dominates its resource area and thus A^ax = '^f'root max^ where the root radius 
is related to maximum size following rroot,max = Crootr^ax Noting that space-filling implies that 
the total area of each size class, A^, is the same for all k, as is the total metabolic rate, , 

requires that A^ = Amax- Combining these relationships gives 

r'k'max 

Atot = / Akdr = T^Croot^max if max ~ ^o) ■ ( 6 ) 

Jro 

which scales as Am oc 


These aggregate quantities allow us to calculate the average biomass density, B: 


D 


Mtot 

Atot 


2ar 

^ ' max 


(7) 


where it is convenient to define ao = l/ctm — 2 q;^. This relationship depends on only two scaling 
exponents, that of the relationship between mass and trunk radius and that of root radius (the lateral 
extent of the roots) to trunk radius. 


In order to make predictions across forests we need to estimate and assess the exponents in 
this relationship which in general can depend on species, environment, forest age, and competitive 
effects (e.g. 14 . 24 - 28 ^ Previous work has shown that theoretically Om ~ 3/8 ^^’^^ando;,. ~ 2/3 

11,15,28,29 consistent with empirical data; together these predict D oc rmax- An alternative way of 
estimating am is to note that many studies report the metabolic scaling relationship with overall 
tree mass, B oc which, together with the relationships described above, gives am = ai/as- 
Theoretical predictions give as — 2 13.30,31 which is in agreement with measurements for a proxy 
of metabolism which give as = 1.78 in one study and as = 1-77 (with a confidence interval of 


6 



1.38 to 2.16) in a second study For the metabolic scaling exponent with mass, ai, there are a 
wide variety of studies and predictions with a common theoretical assessment being that ai ~ 3/4 
(e.g. Ref. Precise measurements of individual tree respiration for a variety of species and 
a diversity of ecosystems (tropical to boreal) have shown a tight relationship between body mass 
and metabolism where the exponent smoothly varied from roughly 1 in small trees to 3/4 in large 
trees, as predicted. Furthermore, given previous analyses of the relationship between tree density 
and average biomass we estimate from Ref. that across modes of plant interaction ai ranges 
between roughly 2/3 and 0.91. Measurements in animals, which we provide as a reference, have 
found that the metabolic exponent varies as widely as 0.25 to 1.25 These results imply that 
could vary in various circumstances along with the resulting scaling for D. For example, taking 
tti between 2/3 and 0.91 and as between 1.78 and 2 predicts a range for of roughly 1/3 to 
1/2 which spans the expected value of 3/8 and agrees with measurements that show a convergence 
to 3/8 as trees become large (with asymptotic behavior at small sizes) Similarly, it has been 
shown that m oc where p is wood density This result would give ~ 3/8 given 

the common assumption that h oc (e.g. ^^) and also shows the importance of species-level 
differences which could affect wood density. Taken together with these calculations predict that 
the biomass density scaling could vary between D oc rliax and D oc rUax- 

To test this result for the dependence of biomass density on maximum size we gathered data 
for the continental United States spanning arid environments to old-growth temperate rain 
forests. Our derivations predict scaling relationships in whole forest properties and so each point 
in our database and plots represent paired values for entire forests where the predicted trends should 
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manifest across forests from many different environments (i.e. many different biomass densities). 
The data represent a variety of forest ages which affects both the distribution of tree sizes and 
how well the observed maximum size represents the true maximum of a given environment. For 
example, previous work shows that forest age, often a measure of recent disturbance, can affect 
both the intercept and exponent of scaling relationships. Using stand age as a proxy for forest 
disturbance we should expect significant error between the predicted scaling relationships and 
observations for young forests suggesting that we examine the scaling of biomass density within 
forest age categories We find that as forests mature the observed scaling exponents between 
biomass density and maximum observed trunk radius converge to the theoretical prediction of 4/3 
with some variation that could be the result of species-, environment-, or forest-level differences in 
the underlying scaling exponents (Figure 1(a) and 1(d), and Table 1 which provides information 
for each age bin and fit). These results illustrate the predictive power of the maximum tree size for 
overall forest structure and also the strong effect that forest maturity has on the deviations between 
observed and steady-state forest properties. In Figure 1(d) we have also shown other possibilities 
for an based on the above calculations. 


Similarly, it has been shown that the net primary productivity of a forest scales with plant 
mass in roughly the same way that metabolism does and thus the two are linearly related: 
NPP = CpB, where Cp is a constant. The average NPP is then given by 


D 

NPP = C OC r‘^B-2ar 

1 1 A ^ • max 

^tot 


( 8 ) 


where we define aNPP = oib — 20;^- The theoretical values discussed earlier for the exponents 


would predict that NPP oc rUax. Again considering possible variation in the exponents, we find 
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NPP oc when we take as = 1.77. 


Data show that these predicted scaling relationships also hold for forests that reach 
a certain age of maturity (Figure 1(b) and 1(e)). Taken together our relationships for biomass 
density and primary productivity imply that they are related via 

NPP (X £)i°‘B-2ar)/{l/a m ) (9) 

which, using predicted exponents, would give NPP oc This is well supported by data from 
a range of forests analyzed here (Figure 1(c) and 1(f)) as well as in previously published data 

(Figure 3). 

It should be noted that all of the scaling relationships for whole forest averages are con¬ 
sistent with the biomass density and NPP of the maximum size tree, which dominates its re¬ 
source area and should be representative of the overall forest. For example, for the largest tree 
^maxl^max OC , which is consistcnt with the whole forest average biomass density. In 

the above derivations we have allowed each scaling relationship to have an independent exponent 
value for the greatest degree of generality. However, it should also be stressed that many of the 
predicted exponents are related to, or derivative of, other exponents via theoretical derivations and 
empirical results. 

It is interesting to note that for both biomass density, D, and NPP the size-distribution scaling 
exponent, is eliminated in the limit of large trees. For completeness we discuss here the values 
of an which might be needed for calculations of Mtot, Btou or Atot- It has been shown empirically 
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that the size distribution scaling ranges from Ri —3.5 to ~ —2 with young forests clustered 
around —3 and mature forests around —2 A detailed study of forests of different age showed a 
convergence to the theoretically predicted value of —2 as forests mature, similar to the perspectives 
presented in this study. These results along with our own illustrate that many scaling predictions 
may hold only in mature forests with low disturbance. Similarly, many other factors such as species 
differences may affect scaling relationships as discussed earlier and thus caution should be applied 
in the application of any scaling relationship to a particular context. For example, it has been 
shown that the plant density exponent can vary from —1.96 in angiosperm communities to —0.78 
in conifer communities 

As already discussed, recent work has shown how tree allometry can be used to derive the 
energy and water budgets of a tree and to predict maximum tree size as a function of local mete¬ 
orological conditions Because maximum tree size determines the overall density of biomass, 
growth, and metabolic rate within a forest, these previous results provide a simple allometric frame¬ 
work for predicting whole forest features from local resource constraints (Figure 2(a)). Notably 
NPP and biomass density can be written as functions of only maximum size which is then a 
complicated, but calculable, function of environmental factors, 

NPP K (10) 

C « (11) 

Tmax = f {P,t,rh,s) , (12) 

illustrating that the previously published dependencies can be greatly simplified (p is precipita- 
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tion, t is temperature, rh is relative humidity, and s is solar radiation). To illustrate this point we 
used the model in Ref. to predict fmax as a function of only precipitation (holding all other cli¬ 
matic variables to global averages) and from r^ax we are able to predict the dependence of NPP on 
precipitation, in good agreement with published data and our own analysis of the United States 
(Figure 2(b)). We provide this analysis to illustrate how allometric frameworks can be used in the 
future to predict forest features from local climate. It is important to note that in these predictions 
it is necessary to consider the full combination of local climatic variables as increases in one fea¬ 
ture, for example temperature, can have opposing effects in different climates because of differing 
impacts on the tree’s energy and water budgets. This is illustrated by previous studies that show 
differing correlations between temperature and biomass density in different forests and par¬ 
tially explains why Michaletz et al. find only strong correlations with precipitation. Indeed, future 
efforts are needed to understand how NPP and biomass density are related to climate combinations 
in connection with many past studies (e.g. 16 , 18 - 22 , 37 , 38 ^ 

Our work shows that the allometric theories are a powerful tool for predicting forest features 
but that caution should be applied in using these theories as many processes, such as disturbance 
effects, can dramatically alter the observed scaling relationships. It should be noted that if the true 
upper bound on tree height can be predicted then disturbance processes can be calculated from 
the theory presented here. For example, NPP will be constant across all size classes and so the 
NPP calculated for the largest tree should represent the observed value even in a disturbed forest 
providing it fully utilises all existing resources. Similarly for biomass density one can use the 
predicted upper bound to calculate Uk and then integrate equation 4 up to the observed maximum 


11 



size to get the true biomass density of the disturbed forest. In this later ease this will lead to 
lower biomass density in reeently disturbed forests beeause young trees have a higher metabolie 
rate per unit mass (ai < 1). These ealeulations eonneet the steady-state work presented here 
to future efforts to quantify disturbanee proeesses and highlight the importanee of predieting the 
upper bound on tree height whieh is the subjeet of ongoing efforts (e.g. ^5, 39 , 40 ^ Furthermore, 
an underlying assumption of this model is that trees seek to be as tall as possible however this 
is dependent on the types of speeies that live in a given environment and the true upper bound 
on height may not be aehievable if elimate has ehanged reeently and speeies have been unable to 
migrate (e.g. This may be an additional souree of error between the predieted upper bound 

on tree height and observed values of NPP or biomass density. 

It should be noted that most of the perspeetives used within this analysis are based on funda¬ 
mental physieal eonstraints or basie biologieal proeesses sueh as meehanieal stability, spaee filling, 
resouree eompetition, and mortality, and these results represent numerous detailed derivations and 
predietions from a large number of researehers. The final produet, however, is eoneeptually simple 
in that total forest net primary produetivity and biomass density seale with the stem radius of the 
largest tree aeeording to exponents of 2/3 and 4/3. Furthermore, our analysis here shows how var¬ 
ious allometrie perspeetives ean be eombined to prediet forest strueture from loeal resourees. The 
key link is that many forest features ean be derived from maximum tree size whieh in turn is eon- 
strained by meteorologieal eonditions. Moving forward these results provide a redueed perspeetive 
on key Earth-system proeesses sueh as biomass produetion rates and ean be used to foreeast both 
eeologieal response to ehanging elimate and also the biogenie elimate feedbaek. In addition, it will 
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be important to integrate these perspectives with the rich history and current efforts of ecological 
and climate models that are able to handle detailed species-level differences, real-time forest dy¬ 
namics and disturbances, the complex processes of soil dynamics, and the interconnection between 
forest processes and transient meteorology 2 , 4 , 7 , 8 , 44^8 

Methods 

Field-measured stem radius was obtained from the Forest Inventory and Analysis (FIA) data 
spanning years from 2003 to 2007 and representing over two million trees. The aboveground 
biomass density map (circa 2003) was derived from the USDA Forest Service while we obtained 
the NPP estimates from a product (MOD17A3) of the Moderate Resolution Imaging Spectrora- 
diometer (MODIS) averaged from 2001 to 2005 For the forest stand ages, the North American 
Carbon Program (NACP) data were used (circa 2006) This MODIS product measures the sur¬ 
face absorption of Photosynthetically Active Radiation (APAR) and converts the APAR to the NPP 
based on the radiation use efficiency parameter, vegetation type, and climatic conditions (e.g. 
Although the data has been validated with Eddy flux data, ecosystem model simulations, and at¬ 
mospheric CO 2 measurements it should be noted that such spacebome information may have 
errors resulting from aerosol/cloud contamination and saturation in dense vegetation (e.g., ^°). 

All gridded data were resampled at a 1-km spatial resolution, where, for stem radius we 
used a 3 km by 3 km moving window to search for the FIA maximum in order to correct for 
the random “fuzzing” that has been intentionally applied by the Forest Service and introduces 
errors up to 1.6 km Climate data used in this study were derived from the Parameter-elevation 
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Regressions on Independent Slopes Model (PRISM) The PRISM provided long-term mean 
annul total precipitation at the 800 m spatial resolution (climatological averages of 1971-2000 
The ecological region map was implemented to classify the data within each age bin (Table 1). 
We used an eco-regional unit, i.e. province (n = 36). According to Cleland et al. each boundary 
embraces distinctive biotic and environmental factors governing ecosystem structure and function: 
provinces represent climatic subzones in relation to geographical location and vegetation type. 

For each of the observed scaling relationships we found exponents by applying both non¬ 
linear least squares fitting of the power law to the raw data and linear fits to the logarithmically 
transformed data. The results from the two methods were similar (Figure 1 and Figure 4). For 
predictions of fmax, and subsequently NPP, from precipitation we solved for the intersection of 
basal metabolism (Qo) with available water (Qp) as described in Kempes et al. The key param¬ 
eter here is 7, the water absorption efficiency, reflecting the local soil and terrain properties. The 
hypothetical 7 of 1/3 from Kempes et al. 2011 was updated using the topographic wetness index 
(TWI) which takes into account both water-flow direction and accumulation. For the continen¬ 
tal United States, we calculated and normalized the TWI with that of a flat area (relative wetness 
index: rTWI= 1 for a flat area). We found that the mean of the rTWI over the continental United 
States is 0.21, which we use for 7. 

Publication Note 

A previous version of this manuscript was originally submitted to a peer reviewed journal on March 
20, 2015, and following referee comments this current version represents the revisions that have 
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been recently submitted to the same journal on May 30, 2015. 
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Figure 1: Nonlinear fits for the sealing of (a) biomass density and (b) net primary produetiv- 
ity based on the maximum tree size as represented by trunk radius, (c) The relationship between 
biomass density and net primary produetivity. (d)-(f) The eorresponding best fit exponents as a 
function of forest age class compared with the theoretical predictions. The data were analyzed 
within 10-year forest age bins for forests ranging from 10 to 350 years old. For the scaling rela¬ 
tionships 500 randomly selected data points are shown for visual clarity. However, in fitting the 
scaling relationships we used all available data (909,524 grid cells) and each age bin contained an 
average of 25,265 data points, ranging between 121 for the 340-350 year age bin to 102,826 for 
the 60-70 year age bin, all in the continental United States (data from^^^^^’^^). The range of B? 
values for power-law fits from all age classes are given and demonstrate a strong correlation with 
a minimum of B? = .64 across all data. The grey dashed lines represent alternative values for the 
exponents as discussed in the text. 
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Figure 2: (a) A schematic of the framework for predicting bulk forest features from local re¬ 

sources. Predictions from this framework (red line) for the dependence of NPP on annual pre¬ 
cipitation for (b) the data from Michaletz et al. and (c) our own analysis for the United States 
(10,000 randomly selected points are shown for visual clarity). Maximum tree size is predicted 
using the framework presented in Kempes et al. and converted to NPP using the results from 
this paper. The prediction matches the data well with = 0.82 for (b) and = 0.46 for (c). 










Figure 3: The scaling relationships and exponents for biomass and net primary productivity for 
the data presented in Michaletz et al. The scaling exponents approach the predicted value of .5 
for mature forests. 
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Figure 4: Linear fits of logarithmically transformed data for the scaling of (a) biomass density 
and (b) net primary productivity (NPP) based on the maximum tree size as represented by trunk 
radius, (c) The relationship between biomass density and net primary productivity, (d)-(f) The 
corresponding best fit exponents as a function of forest age class compared with the theoretical 
predictions. The data were analyzed within 10 year forest age bins for forests ranging from 10 to 
350 years old. For the scaling relationships 500 randomly selected data points are shown for visual 
clarity. However, in fitting the scaling relationships we used all available data (909,524 grid cells) 
and each age bin contained an average of 25,264 data points, ranging between 121 for the 340-350 
year age bin to 102,826 for the 60-70 year age bin, all from the continental United States (data 
from 32-34,36^ gj-gy dashed lines represent alternative values for the exponents as discussed in 
the text. 




































